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Abstract 


We investigate multivariate integration for a space of infinitely times differen¬ 
tiable functions Ts.u ■= {/ £ C“[0,1]® | ||/||jr,_„ < oo}, where ||/||jr,,„ := 

SUp„=(ai,...,a06Ng IIII/ flLl '= - 




-f and u = {uj}j>i is 


a sequence of positive decreasing weights. Let e(n, s) be the minimal worst-case 
error of all algorithms that use n function values in the s-variate case. We prove 
that for any u and s considered e{n,s) < C(s) exp(—c(s)(log n)^) holds for all 
n, where C{s) and c(s) are constants which may depend on s. Further we show 
that if the weights u decay sufficiently fast then there exist some 1 < p < 2 
and absolute constants C and c such that e{n,s) < C'exp(—c(logn)P) holds for 
all s and n. These bounds are attained by quasi-Monte Carlo integration using 
digital nets. These convergence and tractability results come from those for the 
Walsh space into which Ts^u is embedded. 


Keywords: numerical integration, tractability, super-polynomial convergence, 
quasi-Monte Carlo, Walsh spaces, digital nets 


1. Introduction 

In this paper we approximate the integral on an s-dimensional unit cube 



f{x) dx 


by a quasi-Monte Carlo (QMC) algorithm which uses n function values of the 
form 


1 

An,s{f) ■■= V -f{u) for u £ [0,1)C 
2=1 
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One classical issue is the optimal rate of convergence with respect to n. Another 
important issue is the dependence on the number of variables s, since s can be 
hundreds or more in computational applications. The latter issue is related to 
the notion of tractability if we require no exponential dependence on s. 

A large number of studies have been devoted to numerical integration on 
the unit cube for various function spaces. One typical case is that functions 
are only finitely many times differentiable, e.g., functions with bounded varia¬ 
tion, periodic functions in the Korobov space and non-periodic functions in the 
Sobolev space, see 14, EEi and the references therein. For these cases, it 


is known that the rate of convergence is 0{n~°‘) for some a > 0 and thus we 
have polynomial convergence. Another interesting case is when the functions 
are smooth, i.e., infinitely times differentiable. Dick gave reproducing kernel 
Hilbert spaces based on Taylor series for which higher order QMC rules achieve 
a convergence of 0{n~°‘) with a > 0 arbitrarily larra. The spaces were later 
generalized in (^ . Further results were proved in BII3 , where it is shown 
that exponential convergence holds for the Korobov space of periodic functions 
whose Fourier coefhcients decay exponentially fast. Exponential convergence 
means that the integration error converges as 0{q^^) for some q € (0,1), p > 0. 
Note that exponential convergence was also shown for Hermite spaces on K® 
with exponentially fast decaying Hermite coefficients Q. 

In this paper we focus on a weighted normed space of non-periodic smooth 
functions 


:=UeC“[0,l]® 


■ = 




sup 




Hi 


n s ct- 

i=i ^3 


< oo 


ai“i 


( 1 ) 




■dxs 


rf- It 


with a sequence of positive weights u = {uj}j>i, where := 
is easy to check that all functions in jFs,u are analytic from Taylor’s theorem. 
For s = 1, it is known that the integration error usiM Gaussian quadrature 
with n points converges factorially, see, for instance, [J, (6.53)]. Our interest 
is the multivariate QMC rules using so-called digital nets. This is motivated 
by the results by Yoshiki 12211 and is closely related to the notion of Walsh 
figure of merit (WAFOM) |ld . first introduced by Matsumoto, Saito 

and Matoba. WAFOM is a criterion for numerical integration using digital 
nets and is computable in a reasonable time. Hence we can search for good 
digital nets with respect to WAFOM by computer, see 12, Hi for numerical 
experiments. We observe that generalized WAFOM works well for the space 
J-s,u, see Remark [nm 

The first purpose of this paper is to show that the integration error of QMC 
rules using digital nets for J^s.u can achieve super-polynomial convergence as 
0(exp(—c(logn)^)) x for all s and u considered. Here the hidden 

constant and c may depend on s. We remark that this convergence behavior 
was first observed in [l3| as the decay of the lowest-WAFOM value and that the 
combination of [1^ and implies the convergence result for J's,(i/ 2 )j>i ■ 

We also consider tractability for Tg.u- Let us briefly recall the notion of 
tractability (see 15, 1^ 12 1 foi' more information). Let n(£, s) be the information 
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complexity, i.e., the minimal number n of function values which are required in 
order to approximate the s-variate integration within e. An integration problem 
is said to be tractable if n{£, s) does not grow exponentially in e nor s. In 
particular, two notions of tractability have been mainly considered: polynomial 
tractability, i.e., n{e,s) < and strong polynomial tractability, i.e., 

n{£,s) < for Ti,r 2 > 0. A common way to obtain tractability is to 

consider weighted function spaces as introduced by Sloan and Wozniakowski 
[l9l | . Weighted spaces here mean that the dependence on the successive variables 
can be moderated by weights. Our weights u play the same role. For tractability 
results for spaces of smooth functions, see also [a|. 

The second purpose of this paper is to give a sufficient condition to achieve 
super-polynomial convergence with strong tractability. We show that if the 
weights u decay sufficiently fast then the integration error of QMC rules us¬ 
ing digital nets for J^a,u can achieve dimension-independent super-polynomial 
convergence as 0(exp(—c(logn)P)), where the hidden constant and c are inde¬ 
pendent of s and n, and 1 < p < 2 is determined from the decay of u. This 
implies n(£,s) < Cexp(c(loge“^)^/P) for some C, c > 0. 

These convergence and tractability results are also shown for the so-called 
Walsh space >Vs,a,h into which is embedded, which implies those for J's,u- 
For yVs,a,b, we show that the rate of convergence is of order exp(0(—(logn)^)) 
for all a considered and the strong tractability result is equivalent to a typical 
decay of a. 

The rest of the paper is organized as follows. In Section [2l we give the 
necessary background including Walsh functions, the Dick weight, definitions of 
our function spaces and embeddings among them. In Section [3l we give precise 
definitions of the notions of convergence and tractability used in this paper. 
In Section 21 we present Theorems 14.11 and 14.21 which are the summary of all 
results in this paper. Necessary and sufficient conditions for convergence and 
tractability are given in Sections [5] and [H respectively. 

2. Preliminaries 

Throughout this paper, we shall use the following notation. Let N be the set 
of positive integers and No := N U {0}. For a positive integer b > 2, let Zb be 
a cyclic group with b elements, which we identify with the set {0, 1}, 

equipped with addition modulo b. The operators 0 and 0 denote the digitwise 
addition and subtraction modulo b, respectively. That is, for k, k' S Nq whose 
6-adic expansions are k = i‘^ib^~^ and k' = i^'ib''~^ with S Zb 

for all z, 0 and 0 are defined as 

OO OO 

k®k' = ^77^^-! and A:0fc' = '^r]'^b^-\ 

2—1 i—1 

where rji = Ki + k' (mod b) and r][ = Ki — k' (mod 6), respectively. In case of 
vectors in Ng, the operators 0 and 0 are applied componentwise. 
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2.1. Walsh functions 

In this subsection, we introduce Walsh functions and Walsh coefficients, 
which are widely used in analyzing the integration error, see [^, Appendix A] 
for general information. We first give the definition of Walsh functions for the 
one-dimensional case and then generalize it to the higher-dimensional case. 

Definition 2.1. Let b > 2 be a positive integer and let uJb = exp(27r\/—!/&)• 
We denote the b-adic expansion o/ fc S No by k = ki + K 26 -I- • • • -I- 

with G Zf,. Then the k-th b-adic Walsh function f,walfe: [0,1) —>■ 

{l,a;f,,... , 0 ;^“^} is defined as 

bwalfc(x) := 

for X G [ 0 , 1 ) whose b-adic expansion is given by x = f,ib~^ H— ■, which 

is unique in the sense that infinitely many of the are different from 6 — 1 . 

Definition 2.2. Let 6, s G N with b > 2. Let x = (xi,...,Xs) G [0,1)® 

and k = {ki,... ,ks) G Nq. The k-th b-adic Walsh function bwalfc: [0,1)® —>■ 

{l,a;{,,... , 0 ;^“^} is defined as 

S 

bwa\k{x) := f,walfc^(xj). 

1=1 


Since we shall always use Walsh functions in a fixed base 6 , we omit the sub¬ 
script and simply write walfc or walfc in this paper. Some important properties of 
Walsh functions, used in this paper, are described below, see [1, Appendix A.2] 
for the proof. 


Proposition 2.3. The following holds true: 

1. For all k G Nq, we have 

/■waiax)d», = (‘ = 

Jq I U otfieTwise. 

2 . For all k,l G Nq, we have 



walfc (x) walj {x)dx 


1 if k = I, 

0 otherwise. 


3. For all k, k' G Ng and x G [0,1)®, we have 


walfc®fc'(a;) = walfc (x) walfc-(x), 
walfcefc'(a;) = walfc (x) walfc-(x). 


4. The system {walfc j k G Ng} is a complete orthonormal system in L^[0, 1)® 
for any positive integer s. 
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We define the Walsh coefficients as follows. 

Definition 2.4. Let k G Nq and f: [0,1)'* —?> C. The k-th Walsh coeffieient of 
f is defined as 

/(fc) := [ fix)walk{x)dx. 

J[o,iy 

The Walsh series of the function / is given by 

/(a=) ^ X! /(fe)walfc(a;) 

fcGNg 

for any / G We note that all functions considered in this paper are 

equal to their Walsh series, see also the next subsection. 


2.2. Function spaces and embeddings 

In this subsection, we introduce the function spaces Fs^u, Wg.a.b and Wg.a.b 
considered in this paper and give embeddings from Fs^u to Ws,a,b- 

The space of smooth functions Fs,u is defined as in ©■ Throughout the 
paper, we always assume that 


Wl > U2 > • • • > 0. 


( 2 ) 


It is shown in j22| for b = 2 and [2l| for the general case that Walsh coeffi¬ 
cients of functions in Fs^u decay sufficiently fast, as given in Theorem l2.6l below. 
To state the theorem, we define the Hamming weight v{k), the generalized Dick 
weight pL{a\k) and the modified Dick weight il{a;k) for k G No, a G and 
fc G Nq. Note that the Dick weight is originally defined for a = 0 in 0. 


Definition 2.5. Let a G R and /c G N with b-adie expansion k = 
with Ki G lib- Define the function h as / i ( k ) = 0 for k = 0 and h{K) = 1 for 
K 0. The Hamming weight v(k) is defined as the number of nonzero digits for 
the b-adic expansion of k, i.e., 


v{k) := ^ h{Kt). 

i=l 


We define the generalized Dick weight pL{a; k) and the modified Dick weight 
fl{a; k) for the 1-dimensional case as 


OO CXD 

fj,{a;k) :=''^^{ia)h{Ki) and fi{a; k) := max(i -|- a, l)h{K,i). 

i=l i=l 

For the s-dimensional case, let a = (oi,..., Os) G R® and k = {ki, ..., kg) G Ng 
and we define fx{a; k) and jl{a; k) as 


S S 

pL{a;k):=''^^p,{aj;kj) and fl{a;k) :='^^Jl{aj-,kj). 

1=1 
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We can now give the decay of Walsh coefficients which appears in 
lary 3.10]. 


21, Corol- 


Theorem 2.6. Put := 2sin(7r/&) and Mf, := 2sin([6/2j7r/&). Assume f £ 
J's.u- Then it follows that 


i=i 

where Di, = 2 for b = 2 and Dh = Mf, + bmb/{b — Mb) otherwise. 

This decay motivates us to define Walsh spaces >Vs,a,b and Ws.a.b of Walsh 
series whose Walsh coefficients are controlled by the generalized (resp. modified) 
Dick weight. Let a = {aj)j>i be a sequence of real-valued weights. Throughout 
the paper, we assume 

Ol < 02 < 03 < • • • , (3) 

which corresponds to ©■ We first define Ws,a,b as 


W.,a,b := |/: [0,1)^^ 

equipped with the norm 

IL 

and >Vs,a,& as 


fix) = ^ /(A;)walfc(a;) and 


< oo 




:= sup 
feSN" 


■■={f£ Ws,a,b 


IW, 


:= sup < oo 


feere 


Note that all Walsh series in Ws.a.b and Ws,a,b converge. Indeed, for all X G 
(—1,1) and a positive integer I, we have 


S I 


E 




E 


i+aj)h(K.j^i) 


ke% 
kj <b‘ Vj 


feeNg j=li=l 
kj<b‘ Vj 

s I 6—1 

=nnE 

j—1 i—1 


s I 

j=i i=i 


where we denote the 6-adic expansion of kj by kj = ^ with G Zf, 

in the first equality, and the right-most product converges for Z —>■ oo if jXl < 1. 
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This is also true for the modified Dick weight with ^(a; k) and i + aj replaced 
by /i(a; fc) and max(z + aj, 1). Hence we have 


s oo 


j—1 i—1 

for all AT < 1, 

(4) 

s oo 

Y = Y[ n(l + (^ - 

j—1 i—1 

for all |Ar| < 1. 

(5) 


Thus all functions in Ws.a.b and Ws,a,b converge. 

We now give embeddings from J^s,u to Ws,a,b- From Theorem l2.6l we have 

S 

f=l 

s 

< llyll^ 5-M(-log6(-D6m-i«3);fc) 

f=l 

for / G J-s,u- Thus we obtain continuous embeddings 

C Ws.-u'.b with ll/llw,,„,,, < (6) 

with Il/Ilw..„.., < (7) 

where u' = (—logj(Df,m^^Uj))j>i and u" = (—logf,(m^^Ujy)j>i. Note that 
all functions in Ts,u are equal to their Walsh expansions, see [Sj, Section 3.3] or 
d, Theorem A.20]. Embedding ([B]) implies that good algorithms for Ws,u>,b are 
also good for J^s,u- Thus we mainly consider VVs,a,h in the following sections. 

The Walsh space Wg.aM is considered instead of yVs,a,b in Section [B] since 
the modified Dick weight does not take negative values and thus is easier to 
treat. Actually, >Vs,a,6 and yVs,a,b are norm equivalent. Indeed, we have 


fi{a; fc) < fl{a; fc) < fi{a; fc) + ^ ^ (1 - (i + a^)) 

j=l i&Mj 


for all fc G Ng, where Nj is defined as A/)- := {i G N | i + aj < 1}. Thus 




Iw. „ 6 - 




IWs 


,a.,h ) 


( 8 ) 


where the empty sum equals 0, which implies the norm-equivalence of Ws,a,b and 
y^s,a,b- This means that we can consider >Vs,a,6 instead of Ws.a.b for convergence 
results. Furthermore, in Section 16.31 where we consider tractability results, we 
shall assume a condition on the weights which implies that the constant factor 
in (jSj) is bounded independently of s. 
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3. Integration 


Let % = Ts,u, Ws, 0,6 or Ws,a,b- We consider multivariate integration 


/(/) = / f{x) dx for all / G "H. 

J[o.iy 

We approximate /(/) by algorithms 

An,sif) = , f{tn)), 

where (j)n,s ■ R" —>■ R is an arbitrary function and U G [0,1)® for i = 1,... ,n. 
The worst-case error of the algorithm An^s in the space 'H is defined by 

sup |/(/)-yl„,,(/)|. 
fan 
ll/llw<l 

Let e{n,s,'H) be the n-th minimal worst-case error, 

e(n,s) = e(n,s,'H) = inf 'H), 

where the infimum is extended over all algorithms using n function values. For 
n = 0, we approximate 1(f) by a real number. Since is symmetric, i.e., f €1-1 
implies —f€li, the zero algorithm is the best for n = 0, and thus we have 
e(0, s, H) = 1. Hence the integration problem is well normalized for all s. 

For e G (0,1), we define the information complexity of integration 

n(e, s) = n(e, s, H) = min{n G N | e{n, s, H) < e} 


as the minimal number of function values needed to obtain an ^-approximation. 

We are interested in the convergence of the minimal worst-case error of the 
form 

e(n, s) < foralls,nGN, (9) 

where C{s) and c(s) are positive real numbers which may depend on s and where 
p > 1. The condition p > 1 implies that this convergence is super-polynomial. 
We note that if ([9]) holds, then for all s G N and e G (0,1) we have 


n(e, s) < 


exp 


/ log( 7 ( 5 ) -blogg" 
V c(s) 



( 10 ) 


where (7f)+ := max(X, 0) for G K. Furthermore, (flUl) for e = (7(s)e“°*^®i^^°®"'^^ 
implies dH]). This means that ([5]) is equivalent to (flUl) . which shows that asymp¬ 
totically n(e,s) increases with order exp((7((loge“^)^/P)) with respect to e. 
However, how does n(s,s) depend on s? This, of course, depends on C(s) 
and c(s) and is the subject of tractability. Tractability means that we control 
the behavior of C(s) and c(s) and rule out the cases for which n(e, s) depends 
exponentially on s. In this paper we consider two convergence behaviors as 

e(n, s) < (7exp(Hs)e“‘^^^°®"^^ foralls,nGN, ( 11 ) 
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and 


e(n, s) < Ce fQj- s, n S N, 


( 12 ) 


for some p > 1 and A,C,c > 0. Applying (fTUl) to (fTTl) and using the inequality 
(A + y)i/p < Ai/p + yi/p for A, y > 0, (dH implies that 

n(e,s) < C'exp(c'(loge“^)^'^^)exp(A's^'^^) for all s S N, e G (0,1) (13) 

for some A', C',c' > 0. Conversely, dT51) implies ((TT|) for some A, C,c > 0, which 
follows from the inequality + V^/p) < (A + V)^/p for A, A > 0. 

Similarly m is equivalent to the fact that 

’^(E) s) < C"exp(c'(loge“^)^/^) for all s G N, e G (0,1) (14) 

for some C, c' > 0, which can be regarded as super-polynomial convergence 
with strong tractability. Although ((T5)) shows that the information complexity 
may depend super-polynomially on s, we will show that (ED is equivalent to 
ED for the Walsh space. 

4. Main results 

In this section, we present the main results of this paper. The following the¬ 
orems give the super-polynomial convergence and tractability results for Ws,o,t) 
and ^ s,iA' 

Theorem 4.1. Consider integration defined over the Walsh space Ws,a,b with 
a weight sequence a satisfying (jS]). Then we have the following. 

1. For fixed a, there exist positive constants ci^s (i = 1)2, 3,4) which may 
depend on s and a such that it holds that 



< C3,sexp(-C4,s(logn)^) 


for all s and n. In particular e(n, s, Ws,a.b) is of order exp(0(—(logn)^)). 

2. For any a considered, there do not exist A,C,c > 0 such that (HID holds 
for p = 2. 

3. Let 1 < p < 2. Then the following are equivalent. 

(a) The sequence a satisfies liminfj_>oo > 0, 

(b) There exist constants A,C,c>0 such that for all s,n gN we have 


e{n,s) < Cexp(As)e-"P°s'^)" 


(c) There exist constants C, c > 0 such that for all s, n G N we have 

e{n,s) < 


9 



Theorem 4.2. Consider integration defined over Ts,u with a weight sequence 
u satisfying Q • Then we have the following. 

1. There exist positive constants ctnd C6,s depending on s and u such that 

e{n,s,Fs,u) < C 5 ,s exp(-C6,s(logn)^) for all s,n e N. 

2. Let 1 < p < 2 be a real number. If the weight sequence u satisfies 

liminfj^.oo > 0, then there exist constants C,c>0 

such that 

e{n, s) < for all s,n gN. 

These theorems follow from Theorems l5.3flS751 and Corollaries 16. 101 and 16. 151 


5. Lower bounds 

We prove the following lower bound on e(n, s, Wa,a,b) similarly to 0, Theo¬ 
rem 1], which treats the Korobov space. 

Lemma 5.1. Let A be a finite subset o/Nq. Then for all n < |.4| we have 

e{n,s,Ws,a,b) > { max 
\k,k-eA 

Proof. Take an arbitrary algorithm An,s{f) = 4'n,s{f{ti),---,f{'tn))- Define 
gi{x) = Cfcwalfc(a;) for cj. G C such that gi{ti) = 0 for alH = l,2,...,n. 

Since we have n homogeneous linear equations and |.4| > n unknowns cj., there 
exists a nonzero vector of such c^’s, and we can normalize the c^’s by assuming 
that 

max Icfel = Cfe* = 1 for some k* G A. 
keA 

Define the function 

g 2 {x) := Cgi{x)walk-{x) = C ^ Cfcwalfcefc. (a:), 

keA 

where C is defined as C := (maxj. &M(o;feefe*))-i. Then we have 

ll52||w.,„,,=Cmax|cfc6'^(“^'=e'^*)| 

kGA 

< C'max&'"(“^'“®'“*) < C max = 1, 

k^A k,k*&A 

where || • ,, is naturally extended to complex-valued Walsh series. 

We now define a real-valued function f{x) := {g 2 {x) + ^(a;))/2, where 
g^{x) := g 2 {x). Note that = || 52 ||w .,„,6 since ^(a; kQk*) = yL{a; k*Q 

k) for all k. The norm of / is bounded by 

ll/llw.,„,b < (II52 ||w.,„,6 + II^IIw,,„,6)/2 = ||52||w,,<,,6 < 1- 
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We note that An^sif) = 4'n,s{0, • ■ •, 0) and /(/) = Cck* = C. Hence, 

ein,s,Ws,a^b) > \Iif)-A„M)\ > |C(0,0)|. 

Further we consider the function —/. We have || — f\\ws,a.,b ^ 1 An^s{—f) = 
(j)n,s{0 ,... ,0). Hence, 

e(n, S, Ws,a,b) > IH-f) - An,s{-f)\ >\C + </>n,s(0, . . . , 0)|. 

Combining these two inequalities, we have 

e{n,s,V\>s,a,b) >max(|C'-^„,s(0,...,0)|,|C' + (;i„,s(0, ...,0)|) > C. 

Since this holds for arbitrary algorithm we conclude that e{n,s) > C, as 
claimed. □ 


For a non-negative integer d, we now define 

As,d = {fc G Ng I kj < b'^ for all j = 1, 2,..., s}. 

The cardinality of the set is clearly If aj > 0 holds for all j, then 


max 

k,k*^As,d 


-1 


= I max 

k, k^As.d 


-1 


_ — J^—Y^‘j = -L(d.(d+l)/2+ajd) 


where we use k Q k* G As,d for all fc, k* G As,d for the first equality. This 
implies the following corollary. 

Corollary 5.2. Let d G N and assume aj > 0 for all j. Then we have 
e(n,s,yVs,a,b) > b-^‘3=A‘^^A+A3+iA)d) ^11 n < 6*'^. 


We prove a lower bound for the worst-case error and a necessary condition 
to achieve (HU for Ws,a,b in the following three theorems with the following 
notation. For a; G K, we define expj,{a;) := P and max(a,a:) := (max(aj, . 

Theorem 5.3. Let a = G M® and a' = max(a,0). Then we 

have 


e(n,s, >Vs,a,b) > expj, 


(logn)^ 

2s(log&)2 



logn 
slog 6 



Proof. Let n G N. Put d = [log n/(s log 6)J -I- 1, so that < n < 6®'^. It 

follows from Corollarv l5.2l and the embedding Ws,a',b C Ws,a,b that 


e{n,s,Ws,a,b) > e{n,s,Ws,a',b) > b 


E|=i(dV2-Ha'-Hl/2)d) 


> exp^ 




logn 
slog 5 


-f 1 


f=i 



> 



(logn)^ 

2s(log5)2 



logn 
slog 5 
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which proves the result. 


□ 


Theorem 5.4. For any a considered, there do not exist A,C,c > 0 such that 
(fTTI) holds for p = 2. 

Proof. We will argue by contradiction. Let a' = max(a, 0). Suppose that (fTTI) 
holds for some A,C,c > 0. Then this assumption and Theorem 15.31 imply that 


( {logny 
1 2s(log&)2 



logn 
slog 6 





) < e{n,s,Ws,a,b) < C'exp(yls)e 


holds for all s and n. Taking the limit as n goes to infinity, we have l/(2s log b) > 
c for all s. This is a contradiction. □ 

Theorem 5.5. Consider integration defined over Wg^a.b under ([3]). Let 1 < 
p < 2 and assume that (ED holds for this p and some A,C,c > 0. Put r := 
(p — \)/{2 — p). Then we have 


Uo 

liminf — > 0. 


Proof. Let a' = max(a, 1) and N = N{s) := — 1. Corollary 15.21 implies 

e{N,s,yVs,a',b) > expt ^ L«sJ^ j > expt,(-3s[<J^), 


where we use a' +1/2 < a'^ +1/2 < 5 [a'^J /2 in the second inequality. Combining 
this with the assumption that ED holds for yVs,a,b, for all s we have 

expf,(-3sLa',j2) < C'exp(As)e-“('°s^("»^ 

By taking the logarithm and using < N{s) and C < exp(C's), we have 

-3s(log&)[a/J^ < ((7 +A)s- c((s - l)[a/J log&)^’, 


and thus 


holds for all s. Taking the limit inferior as s goes to infinity, we have 

p_lX l/(2-p) 

^ ^ I 

j->oo J'' ~ 

which implies the desired result. 


> 0 , 


□ 
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6. Upper bounds 


In this section, motivated by 13| and its generalization [20|, we prove the 
existence of good QMC algorithms which achieve super-polynomial convergence 
and tractability in Sections 16.21 and [6.31 respectively. Such QMC algorithms are 
given by digital nets. Digital nets are point sets which have the structure of a 
Zf,-module introduced by Niederreiter, see for instance M- We introduce the 
notion of digital nets in the following subsection. 


6.1. Digital nets 

For a positive integer m and a non-negative integer k with its 6-adic expan¬ 
sion k = we define the m-digit truncated vector tTrn{k) G Z™ as 

trm(fc) = {ki,K2,...,k^)^. 

Definition 6.1. Let Gi,..., Gg G be I x d matrices over Z;, with d < 1. 

Let 0 < k < b'^. For 1 < j < s and 1 < i < I, define yi^kj G Z{, os 

iVl,k,j 5 ■ • • 5 yi,k,j ) — Gjtld{k '), 


where the matrix vector multiplication is over Zf,. Then we define 


— 


yi,k,j 

b 


2/2,fey yi,k,j 

62 F 


G [0,1) 


for I < j < s. In this way we obtain the k-th point Xk = ixk,i, ■ •. ,Xk,s)- VVe 
define P = P(Gi,...,Gs) := {xq, ... ,Xi,d_i} (P is considered as a multiset) 
and call it a digital net over hb with precision I, or simply a digital net. 


The dual net of a digital net, which is defined as follows, plays an important 
role in the subsequent analysis. 

Definition 6.2. Let P = P{Gi, ..., Gg) be a digital net over Zb with precision 
1. The dual net of P, denoted by P^ = P^{Gi,.. . ,Gs), is defined as 

:= {fe = (fci,..., fcg) G Ng I Gjtniki) + • • • + Gjtn(fcg) = 0}. 

The next lemma, which is a slight generalization of [^, Lemma 4.75] to our 
context, connects a digital net with Walsh functions. 

Lemma 6.3. Let P be a digital net over 1b and P^ its dual net. Then we have 


\P\ ^ ^ walfc(a;) 

xGP 


1 ifkGP-^, 
0 otherwise. 


From now on, we consider integration defined over Wg,o,b. We use QMC 
algorithms over digital nets. That is, for a digital net P, we use P{f) := 
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\P\-^E.^P f{x), where we identify the digital net P and the QMC algorithm 
on P. The (signed) integration error of / S yVs,a,b by P is calculated as 

l-Pr^ = 1^1”^ /(fe)walfc(a:) - /(/) 

x^P xGP k^fiQ 

fceNg xeP 

= E - /(o) 

= E ?(*’)■ 

fcG-P-L\{0} 

Hence we have 

|P|-iE/(®)-^(/) < E ^ E ( 15 ) 

xeP fceP-L\{o} ’ keP^\{o} 

Remark 6.4. WAFOM, a criterion for digital nets, is defined as a truncated 
version of th e sum on the rightmost side of (1151) for a — 0 in fPi . [^/ and for 
a = 1 in 12^. [^/. Note that fl{a] k) = fi{a] k) in these cases. Thus the sum (and 
the sum with fi^a; k) replaced by fi{a] k) ) can he regarded as an untruncated 
version of WAFOM generalized by weights a, which has not been considered as 
far as the author knows. ITe can say that Fs,u, TVs.a.b and Ws,a,& are function 
spaces for which WAFOM is a suitable guality criterion. 

We now define the minimal weight of P^ by 

dpi := inf ftia: k). 

fcGP-L\{0} 

Then the rightmost side of (fT^ is bounded by ||/||^ ^ where the 

sum is extended over all fc £ Nq with fl{a; k) > dpi. This argument implies the 
following lemma. 

Lemma 6.5. Let P be a digital net. Then we have 

e'^°^{P,Ws,a,b)< E ( 16 ) 

fcGNg 

/i(o;fc)>5pj_ 


The right-hand side of (TTHl) will be evaluated in the following sections. 

We now prove a lemma which gives the existence of digital nets whose mini¬ 
mal weight is large, which generalizes 13, Proposition 2] and [i^, Proposition 4], 
First we define 

vob,a(M) := |{fc G Ng I /i(a;fe) < M}\. 
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Lemma 6.6. Let M be a real number and pf, be the smallest prime factor of 
b. Let d,l be positive integers with I > M — ai — 1. Lf vo\s.a{M) < Pb‘^ holds, 
then there exists a digital net P over with precision I satisfying |P| = b‘^ and 
dp± > M. 

Proof Let Gi,.. .Gs € be matrices. Recall 

^ G7trz(fci) + --- + GjtrKfcs) =0. 

Thus, for a given 0 fc G with kj < U for all j, we have 

|{(Gi,...,G,) e I k e P^(Gi,..., Gs)}\ < 


where ps is the smallest prime factor of b. Hence we have 

|{(Gi,...,G«) G I min fl{a;k)>M}\ 

fceP-L(Gi,...,G.)\{0} 
kj <b’^ Vj 

= 6*^'-|{(Gi,...,G,) G I min p{a-,k)<M}\ 

k&P^{Gu...,Gs)\{0} 
kj<b‘ Vi 

>bsdi- |{(Gi,...,G«)g (Z;,"")« |fcGP^(Gi,...,G«)}| 

fc#0, kj<b‘ Vi 
p.{a;k)<M 

Thus, if vols,a(AL) < pb‘^ holds, there exists (Gi,..., Gs) G with 

inf /i(a; k) > M. (17) 

fceP-^(Gi,....G.)\{0} 
kj<b‘ Vi 

Furthermore, from the assumption I > M — ai — 1 we have 

min{/x(a; fe) | fc G kj > b^ 3j} > max(l, oi + ^ + 1) > M. (18) 
Combining (HZl) and m, we obtain the result. □ 


6.2. Super-polynomial convergenee results 

In this subsection, we prove super-polynomial convergence for Ps,u, hVs,a,b 
and Ws,a,b- Taking dH]) and ([5]) into account, we have only to consider VVs,a,b- 


First we prove a bound on vols,o(M) along [ll|. Exercise 3(b), p.332] and its 
modifications [l^illSl; which treat the case of a = 0. Since vols,a(Tf) < 1 holds 
if M < 1, we assume that M > 1. We have 


vols,a(M) = ^ 


1 < 


fcGN“ 

ji{cL\k)<M 


E 

k^ni 

fl{a-,k)<M 


j^jl{a-,k) — 


M 


< E fi{a-,k) — M 

kem 
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for all X € (0,1), and the right-most expression is equal to 0^=1 n^i(l + (^ ~ 
l)^max(i-i-a 3 ,i))/^M g]) , gy taking the logarithm on both sides and using 

the well-known inequality log(l -I- X) < X, for all X G (0,1) we have 

S CXD 

logvob,a(M) < ^^(6--bMlogX-^ (19) 

i=i i=l 

We proceed to bound ifa^. > 0, it is equal to X“^+V(l- 

X). Otherwise, we have 

CXD 

E 

2=1 2! 2+aj<l i: 

CXD 

^ E i + E^*' 

2 : i-\-aj <1 2 ' = 1 

= nj+X/{l-X), 

where nj := \J\fj \ = |{i G N | i -I- Oj < 1}|. Thus, in both cases, we obtain 

CXD ^ 

-min(X“0 1). 

2^1 

Applying this inequality to (fT5|l . we have 

logvols,a(M) < (6- l)y^ (uj + -E_ niin(X“0 l)'j -|-MlogX"^ 
j=i V “ / 

S 

< {b-l)'^{nj + (logX-i)-^min(X“^l)) -bMlogX-^ 

1=1 

(20) 

Putting X = 1/exp(y (b — l)s/M) and using min(X“i, 1) < 1, we obtain 
logvols,a(M) < Ns + 2^J{b- l)sM, 

where we define Ns := {b — 1) ’ll- We have thus proved the following. 

Lemma 6.7. For all M >0 we have 

vols^o(M) < exp (^Ns + 2y^{b- l)sM^ . 

We note that Lemma 1^771 and the fact that vols,a(M) < 1 if M < 1 implies 

vob,a(M) < exp (^(^Ns + 2y(&- l)s^ . (21) 


^max(2+aj ,1) _ 


E 




E 
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Now we give a bound on the right-hand side of m- From Lemma O we 
have 

E 

fceNo 2^0 feeNg 

/i(a;fc)>M M+ 2 </i(a;fc) <M+ 2+1 

oo 

< vob,a(M + i + 

i^O 

oo 

< exp (^Ns + 2V(6 - l)s(A^+ * + !)) (22) 

2^0 

for all M > 0. We can easily check ^/x < xj{2\/B) + '/Bj‘2 for all x,B > 0. 
Applying this inequality with x = M + i + 1, the right-hand side of (l22ll is 
bounded by 

OO 

Y exp (^Ns + ^{b-l)s/B{M + 1 + 1) + v^(&- 1 )sb) 

2=^0 

OO 

= 6exp (^Ns + yjih— l)si3^ Y^ exp (6 - l)s/i? - log 6^ {M + i + 1)^ . 

i=0 

Taking B as \/{b — l)s/B = (log5)/2, we obtain a bound on the right-hand 
side of (|T6|1 by 

Y < C',exp(-(log6)M/2), 

jl{a;k)'>M 

where the positive constant Cg is defined by 

Cs = exp (TVs + (log6)/2 + 2(&- l)s/log6) (1 - exp(-(log&)/2))“b 

Hence Lemma 16.51 implies the following lemma. 

Lemma 6.8. Let P be a digital net. Then we have 

e^°''{P,yVs,a,b) < Csexp{-5p±{\ogb)/2). 

Put Cg := Ng + 2y/{b — l)s. From (1^ . the condition of Lemma 16.61 is satis¬ 
fied if exp{Cg'/M) < Ph, which is equivalent to M < (dlog pb/Cg)^. Therefore 
the following bound on the worst-case error follows from Lemmas 16.61 and 16.71 

Theorem 6.9. Let d be a positive integer. Then there exists a digital net P 
over Zf, with precision I with |P| = b’^ and I > {logpb/Cg)^d^ — 1 — ai such that 

gWor(p^ Ws.a,b) < Cg exp {-C’Jd^) , (23) 

where CJ := {log pb)^ {log b)/{2C'g^). 




E 


Ij-lM+i) 
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In particular, e(6'^, s, >Vs,a,b) is bounded by the right-hand side of ( 1 ^ . 
Therefore we have e{n, s,yVs,a,b) < C's exp (—(^"[logbnj^). Thus embeddings 
([6|) and ([8]) imply the following convergence result. 

Corollary 6.10. There exist constants Ci^s and C[ ^ (i = 1,2,3J which depend 
on s and the weights u or a such that for all n we have 

e{n,s,Ws,a.j,) < Ci,s exp (-C( ^.(logn)^) , 
e{n,s,Ws,a,b) < C2,s exp (-C2_s(logn)^) , 
e{n,s,Ts,u) < C3,s exp (-C3 ,,(logn)^) . 


6.3. Tractability results 

We have proved super-polynomial convergence for the function spaces, but 
this convergence depends heavily on s. In this subsection, we prove a tractability 
result under the assumption of the sufficient condition from Theorem 14.II That 
is, let r > 0 and assume that the sequence a satisfies liminfj_,.oo > 0. This 
implies that there exist a positive real number a and a non-negative integer A 
such that 

aj > af for all j > A. (24) 

Hence hereafter we assume (IMl) . Under this assumption, JVj = {i € N j i + 
a-j < 1} is empty for sufficiently large j. Hence the constant factor in (|S]) is 
independent of s and thus the tractability result for Ws,a,6 follows from that 
for Ws.a.b- We also note that (b — 1) is finite and we denote it by N. 

We also show Item [2] of Theorem 14.21 using the embedding ([6]). The following 
arguments are parallel to those in Section 16.21 

First we prove a bound on vols_a(Af) under the assumption (IMl) . We need 
the following lemma to bound X]j=i 

Lemma 6.11. For all 0 < X < 1, we have 

S 

<r-^T{l/r){a\ogX-^)-^/\ 

where r(z) := exp{—t) dt is the Gamma function. 

Proof. Since is a monotonically decreasing function of x, we have 

S pS nOC 

/ X‘^^^dx< / exp(-aa:’'logX-i)(ia;. 

^ JO Jo 


Now we consider the substitution aa;''logX ^ = z. Then we have dx = 
r-i(alogX-i)-i/’-z(i-’'i/’'(iz, and thus 

poo pOO 

/ exp(-aa;’'logX-i)da; = r-i(alogX-i)-i/’' / exp(-z) dz 

Jo Jo 


= r ^r(l/r)(ologX 


U-l/r 


which proves the lemma. 


□ 


Combining (1^01) and Lemma [6.Ill for all X G (0,1) we have 


logvols,a(AL) <{b-l) [ [ rtj + 


U'=i 


logX- 


E 

j=A+l 




logX-1 


M log X 


-1 


< (b-l) 


A 


ir(l/r)c 


JogX-i (logX-i)i+i/ 

Putting X = 1/exp(M“’'/(^'’+^)) and using M > 1, we obtain 
logvob,a(M) < 


+ iV + M log X" 


where ci = (b — 1){A + r ^r(l/r)a + 1. We have thus proved the 

following lemma. 

Lemma 6.12. Assume (1241) . Then for all M > 0 we have 
vob,a(M) < exp _ 

Note that the bound on vols,a(Af) from this lemma is weaker than Lemma lGTTl 
with respect to M but independent of s instead. 

In the following, we bound the right-hand side of (HU) along the lines of 
Section [nm For M > 0 we have 

OO 

< y]exp(ci(M-f i-f (25) 

fL{a\h)'>M 


We can easily check the inequality 


{r+i)/( 2 v+i) < for all a;, B > 0. 

- 2r + \ 2r + l ’ “ 

Applying this inequality with x = M + i + 1, the right-hand side of (1^51) is 
bounded by 


y^exp 


i=0 


Cl 


r -I- 1 
2r + l 


= &exp I Cl 


2r + l 


B-AM -f i -f 1) -f Cl ^ 

2r -I- 1 / 

exp _ log 6^ (M + * -f 1)^ . 


Now we choose B such that 


Cl 


r -I- 1 
2r + l 


B-'' 


log& 


Thus we have a bound on the right-hand side of dT51) as 


^ 5 -A(a;fc) < c2exp(-M(log5)/2), 

fcGNg 

/i(a;fc) > Af 
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where the positive constant C 2 is defined as 


C 2 = exp 


log& Cir / 2ci(r + l) 


2 r- + 1 \ ( 2 r + 1 ) log 6 


1 - exp(-(log 6 )/ 2 )’ 


Hence Lemma 16.51 implies the following lemma. 

Lemma 6.13. Assume (|24p . If P is a digital net, we have 


e^°’^{P,yVs,a,b) < C2exp{-dp±{\ogb)/2). 

Now we prove the existence of good digital nets. By Lemma [ 6 .121 the condi¬ 
tion of Lemma l 6 . 6 l is satisfied if which is equivalent 

to M < (dlogTherefore we have the following bound on the 
worst-case error independent of s. 

Theorem 6.14. Let d € N and put C 3 = ((logpf,)/ci)^^’'+^)/^’’''"^^. Assume 
dm). Then there exists a digital net P over Zf, with precision I with |P| = b‘^ 
and I > — 1 — ai such that 

gWor(p^ < C 2 exp |^_£lj|i^rf(2.+ l)/(r+l)^ (26) 

In particular, s) is bounded by the right-hand side of (l26ll . Therefore we 
have e{n,s,Ws,a,b) < C 2 exp (-C 4 [log^ nj , where C 4 := C 3 (log 6 )/ 2 . 

Thus embeddings dH) and ([H]) imply the following tractability result. 

Corollary 6.15. Assume d24ll . Then there exist constants Ci and C'^ (i = 1,2^ 
which are independent of s such that for all s and n we have 

e{n,s,Ws,a,b) < Cl exp (-Cd(logn)( 2 ’'+i)/C+i)) , 

e{n,s,Ws,a,b) < C2exp (-C''(logn)(2’-+i)/C+i)) . 

Assume that the weight sequence u satisfies liminfj^oo log(u^- > 0. Then 

there exists constants C 3 ,C'^ independent of s such that for all s and n we have 

e{n,s,Ps,^) < C3exp(-C''(logn)(2’'+i)/C+i)) . 
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